Density-functional theory study of the catalytic oxidation of CO over transition metal 
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In recent years due to improvements in calculation methods and increased computer power, it has 
become possible to perform first-principles investigations for "simple" chemical reactions at surfaces. 
We have carried out such studies for the catalytic oxidation of CO at transition metal surfaces, in 
particular, at the ruthenium surface for which unusual behavior compared to other transition metal 
catalysts has been reported. Ifigh gas pressure catalytic reactor experiments have revealed that 
the reaction rate over Ru for oxidizing conditions is the highest of the transition metals considered 
- in contrast, under ultra high vacuum conditions, the rate is by far the lowest. We find that 
important for understanding the pressure dependence of the reaction is the fact that Ru (0001) can 
support high concentrations of oxygen at the surface. Under these conditions, the 0-metal bond is 
atypically weak compared to that at lower coverages. We have investigated a number of possible 
reaction pathways for CO oxidation for the conditions of high oxygen coverages, including scattering 
reactions of gas-phase CO at the oxygen covered surface (Eley-Rideal mechanism) as well as the 
Langmuir-Hinshelwood mechanism involving reaction between adsorbed CO molecules and O atoms. 

PACSnumbers: 82.65.J, 71.10, 82.65.M 
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I. INTRODUCTION 

Carbon monoxide oxidation is one of the most exten- 
sively studied heterogeneous catalytic reactions. This is 
due to both its technological importance (e.g., in car ex- 
haust catalytic converters where the active components 
are transition jnetals such as Pt, Pd, and Rh) and its 
"simplicity"ErLl. A microscopic understanding, however, 
of this most fundamental catalytic reaction is still lack- 
ing. We note, however, that steps in this direction 
recently been made via first-principles calculation^ 
Experimentally, it is difficult to probe the state of the 
reactants during reaction - traditionally, information is 
only available before or after the event. Advances in 
surface science techniques, however, afford new informa- 
tion concerning the behavior of the reactants during the 
reaction process. For example, time-resolved scanning 
tunneling microscopy (STM)lI, tirae-resolved electron 
energy-loss spectroscopy (IREELS)l3, time-resolved in- 
frared spectroscopy (TRIS)lI, and ("fast"[iHgh resolution) 
x-ray photoelectron spectroscopy (XPS)ll3. On aJarger 
scale are low energy electron microscopy (LEEM)li2l and 
photoemission electron microscopy (PEEM)t3. Theoret- 
ically, a significant hindrance has been computational 
limitations, but also, for a realistic description of cer- 
tain reactions, new theoretical developments had to be 
awaited; for example, the generalized gradient approx- 
imation (GGA) for the exchange-correlation functional 
has been shown to be crucial for obtaining accurate ac- 
tivation barriers for hy|drpffen dissociation at metal and 
semiconductor surfacesli3o. 

Past studies performed using surface science tech- 
niques, i.e., on well characterized single-crystal metal sur- 
faces under ultra high vacuum (UHV) conditions, have 
shown that CO oxidation proceeds via the Langmuir- 
Hinshelwood (L-H) mechanism in ■adiich reaction takes 
place between chemisorbed reagentsEI. Recent high gas- 
pressures catalytic reactor experiments, which afford the 
study of chemical reactions under "realistic" high pres- 



sure and temperature conditions, support the assignment 
of the L-H mechanism for Pt, Pd, Rh, and IrEl. For Ru, 
however, somewhat anomalous behavior was found which 
indicated that reaction via scattering of CO molecules 
with adsorbed O atoms may be-|taking place, i.e. via 
an Eley-Rideal- type mechanismll3. In particular, with 
pressures of about 10 torr and for oxidizing conditions 
(i.e., at CO/O2 pressure ratios < 1) the rate of CO2 
production was found to be significaujtljj higher than 
at the other transition metal surfacesli3ii3; in contrast, 
under UHV jjopditions, the rate is extremely low over 
Ru (OOOl)y'Ejy. Unlike the other transition metals, al- 
most no chemisorbed CO could be detected either dur- 
ing or after the reaction, and the kinetic data (activa- 
tion energy and pressure dependencies) was found to be 
markedly different to that of the other metals; in particu- 
lar, highest rates occurred for high concentrations of oxy- 
gen at the surface, whereas for the other metals, highest 
rates occurred for low O coverages. Our studies show that 
Ru does behave differently to the other transition metal 
catalysts in that high coverages of O can be supported 
on the surface (up to desorption temperatures) where the 
0-metal bond is significantly weaker than in the lower 
coverage phases. Investigation of the energetics for CO2 
formation indicates that a Langmuir-Hinshelwood mech- 
anism, rather than an Eley-Rideal process is dominant. 



II. CALCULATION METHOD 

In order gain understanding into the apparently dif- 
ferent behavior of Ru for the CO oxidation reaction, 
and to obtain a microscopic picture of this basic sur- 
face catalyzed reaction in general, we carried out density- 
functional theory (DFT) calculationsEfl. We use the ab 
initio pseudopotential plane wave method and. the su- 
percell approach where we employ the GGAl3 for the 
exchange-correlation interaction. We use ab initio, fnHj 
separable, norm-conserving GGA pscudopotentialsE^" 



where for the Ru atoms, relativistic effects are taken 
into account using weighted spin-averaged pseudopoten- 
tials. The surface is modehed using a (2 x 2) surface unit 
ceU with four layers of Ru(OOOl). An energy cut-off of 
40 Ry is taken with threcj-apecial k-points in the two- 
dimensional Brillouin zonec3. The_adsorbate structures 
are created on one side of the slabtj. We relax the posi- 
tion of the atoms, keeping the Ru atoms in the bottom 
two layers fixed at their bulk-like positions. 

III. OXYGEN ON RUTHENIUM 

Under UHV conditions, at room temperature, dissocia- 
tive adsorption of O2 results in an (apparent) saturation 
coverage of Go ~ 1/2, cprresponding to the formation 
of a (2 X 1)=D structurekil. At 60=1/4, a (2 x 2)-0 
phase formsEj. Here Qo = 1 means that there are as 
many O atoms as there are Ru atoms in the top layer. 
In both surface structures, O adsorbs in the hexagonal 
close-packed (hep) site. p.ux. earlier DFT-GGA calcula- 
tions for O on Ru (OOOl)BEl indicated that even higher 
coverage phases should form; namely, a (1 x l)-0 struc- 
ture with coverage Go = 1, as well as a (2 x 2)-30 struc- 
ture with coverage Go = 3/4. As for the lower coverage 
structures, the O atoms occupy hep sites. The adsorption 
energy of O decreases notably with increasing coverage 
and for the monolayer coverage phase, the adsorption 
energy is wO.7 eV less than that of the (2 x 2)-0 phase. 
Both the (2 x 2)-30 and (1 x l)-0 structures have subse- 
quently been created experimejatally with the use of NO2 
or high gas pressures of 02EII~e3, and the atomic struc- 
ture verified by dynamical low-energy electron diffraction 
(LEED) intensity analyses. NO2 readily dissociates at el- 
evated temperatures in the presence of adsorbed oxygen 
delivering atomic oxygen to the surface while NO des- 
orbs. Furthermore, it has been demonstrated that after 
completion of the monolayer oxygen strU| Ct\ir e, additional 
oxygen can enter the subsurface regioncZlO at elevated 
temperatures. Formation of the higher coverage phases 
(G =3/4 and 1.0) from gas-phase O2 at "usual" expo- 
sures under UHV conditions is apparently kinetically hin- 
dered by activation barriers for O2 dissociation, induced 
by the pre-adsorbed oxygen atoms at coverage Go ~0.5. 
With respect to the high pressure catalytic reactor exper- 
iments mentioned above, because the conditions under 
which the highest rates of CO2 formation were reported 
involved high O2 partial gas pressures and oxidizing con- 
ditions, there will be a significant attempt frequency of 
O2 to overcome activation barriers for dissociative ad- 
sorption. Thus, it is likely that during reaction, the oxy- 
gen coverage on the surface approaches one monolayer. 
We note also that in the catalytic reactor experiments it 
may be unlikely that there is significant subsurface oxy- 
gen present since the temperature range studied in the 
experiment of 380 K to 500 K is less than that of 600 K 
at which oxygen is reported to enter the. S].ihsurface re- 
gion with an appreciable rate (using N02)EZlO and auger 
electron spectroscopy indicated a coverage of about one 



monolayer. 

IV. REACTION VIA GAS-PHASE CO WITH 
ADSORBED O 

In earlier publications we reported our investigation feu: 
reaction via scattering of gas-phase CO with adsorbed Ocl 
so here we only briefly describe the results. 

The (1 X l)-0 phase is assumed to cover the whole 
surface and we investigate the interaction of CO with this 
oxygen-covered surface. For a given lateral position, CO 
is placed well above the surface (with the C-end down) 
and the total energy calculated for decreasing distances 
of CO from the surface. All atomic positions are relaxed 
except that of the C atom which is held fixed (and the 
bottom two Ru layers). 




FIG. 1. (a) Schematic diagram illustrating the on-top (t), 
fee (f), hep (h), and bridge (b) sites of (1 x l)-0/Ru (0001). 
Small and large eireles represent O and Ru atoms, respec- 
tively, (b) Top view of CO adsorbed in an 0-vacancy of 
(1 X l)-0/Ru (0001). The full lines represent the (symmetry 
inequivalent) region within which we consider reaction takes 
place. CO is indicated by the small black circle. 

We considered a number of lateral positions: the on- 
top and fee sites, with respect to the Ru (0001) substrate, 
a bridge site between two adsorbed O atoms, as well as 
directly above an adsorbed 0(a) atom (see Fig. la). We 
find that an energy barrier begins to build up at about 
2.5 A from the surface (with respect to the average po- 
sition of the O atoms) for all sites, reflecting a repulsive 
interaction of CO with the 0-covered surface. A very 
weak physisorption well of about 0.04 eV is also found 
above the 0-adlayer. Its position is Ri3.0 A from the 
surface for the on-top and fee sites, and at about !v3.5 A 
for the bridge site and the approach directly over the O 
atom. Thus, for a full monolayer coverage of oxygen on 
the surface, CO is unable to form a chemical bond with 
the substrate and the L-H reaction mechanism is there- 
fore prevented. 

In order to obtain a more detailed understanding of 
CO2 formation via a scattering reaction, we evaluated an 
appropriate cut through the high-dimensional potential 
energy surface (PES) (see Ref.a). This cut is defined 
by two variables: the vertical position of the C atom and 
the vertical position of the reacting 0(a) adatom directly 
below. Initially, with the CO-axis is held perpendicular 
to the surface we find that the activation barrier for CO2 
formation is «1.6 eV. When the tilt angle of the CO- 
axis is allowed to relax, the energy barrier is reduced to 
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about 1.1 eV. The transition state (depicted in the inset 
of Fig. 2) has a 
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FIG. 2. Calculated energy diagram for the E-R mechanism 
of CO oxidation at Ru(OOOl). Note that the depths of the 
physisorption wells are exaggerated for clarity. The transition 
state geometry is indicated in the inset. The large, medium, 
and small (dark) circles represent Ru, O, and C atoms, re- 
spectively. 

"bond angle" of « 131° and the reacting 0(a) atom 
is «0.35 A above the other O atoms in the surface unit 
cell. The C-O(a) bond length is 1.50 A (stretched by 
27 % compared to the calculated bond length of a free 
CO2 rnolecule which is 1.18 A ; the experimental value is 
1.16 AeI), and the bond length of CO is 1.17 A (the cal- 
culated value of a free CO niolecule is 1.15 A ; the exper- 
imental value is also 1.15 Al3). After onset of reaction, 
the CO-O(a) bond begins to develop and the 0(a)-Ru 
bond is weakened. It then becomes energetically unfa- 
vorable for the reacting complex to be at the surface and 
it is strongly repelled towards the vacuum region. The 
resulting energy diagram is shown in Fig. 2. It can be 
seen that there is a significant energy gain from the sur- 
face reaction of 1.95 eV. 

Using the determined activation barrier in a sim- 
ple Arrhenius-type equation with the prefactor obtained 
from consideration of the number of CO molecules hitting 
the surface per site per second at a given temperature and 
pressure, we can estimate the reaction ratea. This will 
give an upper bound since for other orientations of the 
molecule the barrier is larger. The rate is found to be 
significantly Jawer (by 3 x 10~^) than that measured ex- 
perimentallyMl. This indicates that this mechanism alone 
cannot explain the enhanced CO2 turnover frequency as 
was speculated. To investigate other possible reaction 
channels, we consider it conceivable that there |-are va- 
cancies in the perfect (1 x l)-0 adlayer (see Ref.cl for an 
estimate of the 0-vacancy density) . CO molecules may 
then adsorb at these vacant sites and react via a L-H 
mechanism. In the following section we investigate such 
a L-H reaction mechanism. 



V. REACTION BETWEEN ADSORBED CO AND 

O 

We first consider the energetics of adsorption of CO 
into a vacant hep site of the monolayer oxygen structure. 



Interestingly, we find that there is an activation barrier 
of f«0.3 eV. We expect, however, that this barrier will rel- 
atively easily be overcome at the high gas pressures used 
in the catalytic reactor experiments. From Fig. lb it can 
be seen that CO is closely surrounded by a hexagonal ar- 
rangement of six O atoms. Even though the reactants are 
already very close, there is still a repulsive rather than an 
attractive interaction between CO and the neighboring O 
atoms. In determining the reaction path we consider re- 
action in the inequivalent area of the surface, indicated in 
Fig. lb by the continuous lines. There are clearly a num- 
ber of possible ways that the reaction between CO and 
a neighboring O atom could proceed: for example, the 
CO molecule may approach an O atom, an O atom may 
approach the CO molecule, or the reactants may both 
move towards each other over one Ru atom. For each of 
these scenarios there are obviously also a number of pos- 
sible reaction paths (i.e., via the on-top site, bridge and 
fee site, etc). To determine the minimum energy path- 
way, energy barrier, and associated transition state, is 
clearly not a simple problem. Recent ab initio studies of 
surface reactions and dissociation of diatomic molecules 
at surfaces have attempted "direct ™«thods" for finding 
the lowest energy reaction pathwaj^S. In the present 
work, however, we use the "standard grid approach" of 
constructing various relevant PESs since we are also in- 
terested in the shape of the PES away from the minimum 
energy reaction pathway. 

In view of the weaker CO-metal bond strength com- 
pared to that of the 0-metal bond strength at this cov- 
erage, i.e. 0.85 eV compared to 2.09 eV (with respect to 
gas phase 1/2 O2), we first consider reaction via move- 
ment of CO towards the O atom. This in fact turns out to 
have the lowest energy pathway of the ones we considered 
with a barrier of w 1.51 eV. This value is consistent with 
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the recent experimental estimate of >1. 
of high oxygen coverages on the surfacec 
reaction pathways considered were for O to move towards 
CO and for O and CO to move towards each other. The 
lowest energy pathway found for these scenarios were at 
least 0.2 and 0.3 eV higher, respectively. We point out 
that in this work we have only considered the most obvi- 
ous reaction paths, and in order to explore in more detail 
the very complex nature of this high-dimensional poten- 
tial energy surface, further calculations are required. 

As a first step we investigate the energetics for differ- 
ent fixed lateral positions of the C atom within the area 
shown in Fig. lb. We initially keep the lateral position of 
the O atoms fixed, but allow vertical relaxations. For CO 
moving towards the on-top site (see Fig. 1), a strong re- 
pulsion between C and the two symmetrically equivalent 
O atoms develops giving rise to a large energy barrier 
of 2.53 eV. Thus, the pathway over the on-top site is 
energetically unfavorable. 

Interestingly, the situation is somewhat different for 
CO moving towards the fee site; in this direction there 
is also the build up of an energy barrier; on overcoming 
the barrier, however, there is an attractive interaction 
between C and the two symmetrically equivalent O atoms 
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(see Fig. lb). When the O atoms are then aUowed to 
laterally relax, we find the formation of a carbonate-like 
species. The atomic geometry is depicted in Fig. 3a. We 
note that in our investigations of the CO-gas 

b(CO)=1.20A b(CO)=1.19A 
b(C-Oad)= 1 .42 A b(C-Oad)= 1 .42 A 

b(0-Ru)=2.79, 2.28, 2.23 A b(0-Ru)=2.62, 2.41, 2.18 A 

FIG. 3. Atc^mic geometry of ide/ftlfied carbonate species at 
(a) a vacancy in the monolayer oxygen structure, and (b) at 
the perfect (1 x l)-0/Ru (0001) surface. The large, small, and 
small dark circles represent Ru, O, and C atoms, respectively. 
The various bond lengths are indicated (in A). 

scattering reaction, we also identified the stability of 
such a carbonate species on the fully 0-covered sur- 
face; the atomic geometry is similar as can be seen from 
Fig. 3b. Formation of this species was also found to in- 
volve a significant energy barrier. We note that experi- 
mental identification of carbonate species in CO oxida- 
tion reactions ovet-other transition metals have been re- 
ported (e.g., Ref.Efl), where they may possibly act as an 
intermediary. The actual role they play in the carbon 
monoxide oxidation reaction for this system, however, is 
at present unknown. 
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FIG. 4. Calculated energy diagram for the L-H mechanism 
of CO oxidation at Ru (0001). The transition state geometry 
is indicated in the inset. The large, small, and small darker 
circles represent Ru, O, and C atoms, respectively. 

We find that the minimum energy pathway found 
for CO2 formation corresponds to one where the CO 
molecule moves essentially directly towards the O atom. 
The associated energetics are shown ja Fig. 4 where we 
have constructed the energy diagramEJ. Similarly to our 
study of the E-R mechanism, we also find a small ph- 
ysisorption well for CO2 above the surface. The cor- 
responding transition state geometry is depicted in the 
inset of Fig. 4. In this geometry, the C-O(a) bond is 
almost parallel to the surface and the CO axis is bent 
away from 0(a) yielding a bent CO-O(a) complex with 
a bond angle of 125°; similar to that found for the E-R 
mechanism. The C-O(a) bond length is 1.59 A (about 
35 % stretched compared to that in CO2) and the CO 



bond length is 1.18 A. 




FIG. 5. Atomic positions a selected points along the reac- 
tion energy pathway CO2 formation. The large, small, and 
small darker circles represent Ru, O, and C atoms, respec- 
tively. 

Figure 5 shows the atomic geometry at selected posi- 
tions along the reaction path to CO2 formation. Initially 
CO is in the vacancy and 0(a) in the hep site. As CO 
approaches 0(a), repulsive forces build up on both the 
C and 0(a) atoms and the molecular axis of CO begins 
to tilt away from 0(a). At the transition state, CO and 
O begin to lift off the surface as they break their metal 
bonds in favor of developing a C-O(a) bond. Interest- 
ingly, the distance that O moves away from the surface 
by is ~0.36 A ; very similar to that found in the E-R 
mechanism which was «0.35 A . As CO2 be gins to form, 
the molecular axis quickly straightens out to its linear 
geometry. 

The corresponding valence electron density and density 
difference distributions are shown in Fig. 6. The latter 
is constructed by subtracting from the electron density 
of the CO,0/Ru (0001) system, those of the 0-covered 
surface and a free CO molecule. From the electron den- 
sity difference distributions, with respect to the starting 
configuration of CO in the vacancy, it can be seen that as 
CO moves towards 0(a), firstly there is less electron den- 
sity in the CO 2tt* orbitals, indicating a weakening of the 
C-metal bond. Also, electron density has been depleted 
from the 0(a) orbital pointing towards the CO molecule 
and an increase occurs into orbitals in the orthogonal di- 
rection (pointing in a near-perpendicular direction to the 
surface). The redistribution is an effect of Pauli repul- 
sion. At the transition state, the onset of bond formation 
can be seen between these latter 0(a) orbitals and the 
27r*-like orbitals of the C atom. We can notice that CO 
is bonded only weakly to the metal through a 27r*-like or- 
bital. The bond of the reacting 0(a) to the surface is also 
significantly weakened at the transition state as can be 
see from the lower panel showing the total valence elec- 
tron density. Interestingly, the density difference plot for 
the transition state is very similar to that of the E-R 
process (see Fig. 7 of RefJj). In the last panel (leftmost), 
significant accumulation of electron density can be seen 
between the C and 0(a) atoms as the CO2 molecule is 
practically formed. 

To summarize, our studies indicate that a Langmuir- 
Hinshclwood mechanism rather than an Eley-Rideal 
mechanism is the dominant reaction process giving rise 
to the reported increase in reactivity of ruthenium for the 
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CO oxidation reaction as measured in the high pressure 
catalytic reactor experiments as compared to under UHV 
conditions. We find that for high coverages of O on the 
surface, which are attainable under the sufficiently high 
oxygen pressures used in the experiment, the O adsorp- 
tion energy is notably weaker than the lower coverage 
phases that form under UHV conditions. Furthermore, 
the adsorption energy of CO in the presence of high O 
coverages is weaker than in the presence of low O cov- 
erages. These factors, together with the close proximity 
of the reactants is thought to give rise to the observed 
increase in reactivity of Ru and the reported anomalous 
behavior on partial gas pressure. 

Finally, we like to mention very recent results of CO 
oxidation experiments: For the case of very high concen- 
trations of oxygen at the surface (one monolayer on the 
surface plus oxygen occupying subsurface sites) which can 
be prepared after formation of the (1x1) phase by using 
either NO2 or high gas pressures of O2 at eteuated temper- 
atures, reaction rates notably greater than that from the 
on-suxfe££ monolayer oxygen structure have been mea- 
suredEj'EZl. These high rates have been proposed to be 
connected to the existence of copious amounts of oxygen 
in the subsurface region. These are clearly very inter- 
esting results and require more detailed investigations in 
order to understand this beha 'ior. 



A. Conclusion 

We have performed density functional theory calcula- 
tions in order to investigate the catalytic oxidation of car- 
bon monoxide over Ru (0001) for the conditions of high 
oxygen coverages on the surface for which highest reac- 
tion rates of CO2 formation have been reported. It is 
only by exposure of the surface to high O2 gas pressures, 
or with the use of strongly oxidative molecules such as 
NO2 that high oxygen coverages can be achieved. In this 
case the 0-metal bond strength is notably weaker com- 
pared to the lower coverage structures that form under 
UHV conditions, and as such, is expected to be more re- 
active. In the present work we concentrated mainly on 
the microscopic description and study of the Langmuir- 
Hinshelwood reaction mechanism. We identified a bent 
transition state for CO movement towards an adsorbed O 
atom with an associated energy barrier of approximately 
1.5 eV. In this configuration the CO and the adsorbed 
O atom have substantially weakened (and significantly 
stretched) their bond to the substrate. 




FIG. 6. Valence electron density (lower panels) and density 
difference distributions (upper panels) for selected positions 
along the reaction path. The units are in e bohr^'^. For the 
upper leftmost panel the contours differ: The first negative 
contour is at —3 x 10"'^ and the spacing is — 15 x 10"'^ and 
the first positive contour is at 3 x 10~^ and the spacing is 
30 X 10"^ 
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